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We investigate Yang-Lee zeros of grand partition functions as truncated fugacity polynomials of 
which coefficients are given by the canonical partition functions Z{T, V, N) up to Y < Ymax- Such a 
partition function can be inevitably obtained from the net-baryon number multiplicity distribution 
in relativistic heavy ion collisions, where the number of the event beyond Ymax has insufficient 
statistics, as well as canonical approaches in lattice QCD. We use a chiral random matrix model as 
a solvable model for chiral phase transition in QCD and show that the closest edge of the distribution 
to real chemical potential axis is stable against cutting the tail of the multiplicity distribution. The 
similar behavior is also found in lattice QCD at finite temperature for Roberge-Weiss transition. 

In contrast, such a stability is found to be absent in the Skellam distribution which does not have 
phase transition. We compare the number of Ymax to obtain the stable Yang-Lee zeros with those 
of critical higher order cumulants. 

PACS numbers: 12.38.Gc, 12.38.Mh, 25.75.Nq, 25.75.Gz 


I. INTRODUCTION 

Phase transition in quantum chromodynamics (QCD) 
is one of the central subjects in high energy nuclear 
physics both theoretically and experimentally. First 
principle lattice QCD (LQCD) calculations have shown 
that the transition from quark-gluon plasma (QGP) to 
hadronic matter is of crossover type at physical quark 
masses [I|, in which order parameters and thermody¬ 
namic quantities change smoothly as functions of tem¬ 
perature. At finite baryon density, one expects that the 
nature of the transition can change. Unfortunately, lit¬ 
tle is known about the state of matter at high baryon 
density from LQCD calculations because of the difficulty 
in numerical simulation at finite baryon chemical poten¬ 
tial /I UM- Various approximation methods applied 
so far seem to work only /r < T or a small volume or 
heavy quark mass region. Nevertheless, effective models 
which implement relevant symmetries in QCD and large 
Nc studies have shown that rich phase structure exists in 
high density ill- In particular, if there is a first order 
phase transition at T = 0 and large /i, there must be a 
critical point (CP) at which the first order phase tran¬ 
sition line terminates and the transition becomes second 
order. Existence of CP is supported by many chiral ef¬ 
fective models Jg , but its location depends on the detail 
of the models 0. 

Stimulated by these theoretical results, the first beam 
energy scan program at Relativistic Heavy Ion Collider 
(RHIC) has been carried out in search for the CP ii. 
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Since lower colliding energies leaves the incident nucle¬ 
ons in the central region, one expects to explore higher 
baryon density region at lower energies. There are a 
number of observables which might have potential to 
indicate the transition from QGP to hadronic matter. 
Among them, event-by-event fluctuations of conserved 
charges are intimately connected to critical behavior as¬ 
sociated with the phase transition [lol - [l3l |. Measure¬ 
ments of the net-proton number fluctuations as a proxy 
of the net-baryon one 0 and net electric charges have 
been presented for Au-fAu collisions at various energies 
from =7.7 GeV to 200 GeV [lM3. Given 

the fact that multiplicity of different particle species are 
well described by statistical models [isl - l^ . one may re¬ 
gard event-by-event fluctuations of conserved charges as 
those of the grand canonical ensembles at chemical freeze- 
out temperature T and baryonic chemical potential fi. 
Through systematic analyses of the location of (T, /r) cor¬ 
responding to each colliding systems [2^, one can map 
experimental measurements for property of the matter 
on T — plane. Furthermore, recent LQCD results at 
physical quark masses indicate that the crossover region 
coincides with the chemical freeze-out at least for fi < T 
. One may look for remnant of the chiral critical¬ 
ity in the crossover region, originating from second order 
phase transition in the vanishing quark mass [23l - l^ . 

Property of the transition can be characterized by be¬ 
havior of fluctuations of conserved c harges as well as an 
order parameter and its fluctuations [261127|| . In the case 
of the chiral phase transition in QCD, the chiral order 
parameter or quark-antiquark condensate (qq) couples 
to quarks carrying the baryon number and the electric 
charge. Thus, the second order chiral phase transition 
in the chiral limit at finite temperature is characterized 
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by not only divergent fluctuation of the order parame¬ 
ter but also higher order cumulants of the net baryon 
number and the net electric charge [ 2 ^. The divergence 
of the conserved charge fluctuations is governed by the 
critical exponent of the specific heat a which depends 
on the universality class QCD belon gs t o. Although it 
is not completely determined yet |29l . l30| , recent simula¬ 
tions [3l| indicate 0(4) in the three dimensions, as con¬ 
jectured by Pisarski and Wilczek [s^. In this case, the 
first divergent cumulant appears at the sixth order. At 
finite but small quark masses, the divergence is replaced 
by sign change, owing to the property of the universal 
0(4) scaling function [23, HI]. At nonzero net baryon 
number density, the divergence at 2n-th order cumulants 
in the chiral limit appears at n-th order one. The tricrit- 
ical point in the chiral limit becomes the CP, where the 
second order cumulants diverge [23,[13,|35|. The chemical 
freeze-out line may locate at lower temperature than the 
chiral phase boundary such that measured fluctua¬ 
tions might not reflect those at the phase transition [s^ ■ 
Nevertheless, the existence of the CP is accompanied by 
anomalous behavior of the cumulants such as negative 
third and fourth order cumulants around the CP |38l - l40l| 
and may lead to non-monotonic behavior of the higher 
order cumulants as functions of ^ysNN- Indeed, the mea¬ 
sured net-proton number cumulants in seems to fol¬ 
low this expectation, although still inconclusive due to 
uncertainty. 

The measurement of the cumulants is based on event- 
by-event multiplicity distribution. Once the fluctuations 
are regarded as those of the grand canonical ensemble, 
the multiplicity distribution can be identihed with un¬ 
normalized probability distribution. 

While the cumulants are expressed by central moments 
of the probability distribution, it is convenient for theo¬ 
retical studies to compute them by differentiating the 
thermodynamic pressure with respect to chemical poten¬ 
tials. Recently, one of the authors (K.M.) investigated 
the probability distribution of the net baryon number in 
models with phase transitions [^ . 1^ . . It turns out 

that sufficient information on the tail in the probability 
distribution is responsible for the critical behavior of the 
higher order cumulants and that the remnant of the 0(4) 
criticality can be characterized by narrower tail than the 
corresponding reference distribution. In the probability 
distribution, such information on the phase transition is 
encoded in the N dependence of the canonical partition 
function Z (T, V, N). 

Since the grand partition function is more straight¬ 
forward in relativistic quantum held theories where the 
number of particles are not dehnite, computations of 
the canonical partition function are not generally easy. 
In Ref. H^, Hasenfratz and Toussaint proposed that 
the canonical partition function, Z (T, V, N ), is calcu¬ 
lated through the Fourier transformation of the grand 
canonical partition function,2^(T, R, /i), evaluated at pure 
imaginary fj,. The difficulty associated with the complex 
fermion determinant is replaced by the highly oscillating 


inte gra l which requires extraordinary numerical precision 

EUlSlIil. 

Nevertheless, the probability distribution gives further 
insights into property of the system including phase tran¬ 
sitions. 

In Ref. [ 4 ^, one of the authors (A.N.) pointed out 
that one can extract the fugacity parameter A = at 
the chemical freeze-out and construct Z (T, V, N) for the 
net baryon number without any assumption on the prop¬ 
erty of equilibrium P{N). Furthermore, once Z{T, V, N) 
is known, one can obtain the grand partition function 
Z (T, V, fi) as a series of fugacity. This enables us to apply 
Yang-Lee theory for the phase transition [dj (I^or re¬ 
cent reviews, see, e.g., in which zeros of the par¬ 

tition function give information on the thermodynamic 
property of the system. The zeros of the partition func¬ 
tion are distributed on a line in the complex plane of an 
external parameter and its density grows up with the sys¬ 
tem volume, then finally coalesce into the line in the ther¬ 
modynamic limit. This property leads us, in principle, 
to obtain the location and order of the phase transition 
from the distribution of the zeros. Even in the absence 
of the phase transition, the zeros accumulated on the 
edge of the distribution exhibit singular behavior. This 
singularity, known as Yang-Lee edge singularity [Slj, , 
can be regarded as a CP in the complex plane and gives 
influence on the thermodynamics on the real axis |53l |. 

In both experiments and the canonical approach in 
LQCD, Z{T,V,N) at large N requires such high statis¬ 
tics that obtained information is limited to some finite N, 
thus one has to truncate the fugacity polynomial there in 
reconstructing the grand partition function (See [H, 
for recent LQCD calculations). It is not a priori clear 
whether one can obtain the correct information on the 
phase transition from such a truncated partition func¬ 
tion. The purpose of this paper is to clarify this point. 
We employ a solvable model for the chiral phase transi¬ 
tion in QCD. We present the Yang-Lee zeros in a chiral 
random matrix model, both for the exact grand partition 
function and for the reconstructed one as a truncated fu¬ 
gacity series with the canonical partition function being 
the coefficients. We discuss effects of the truncation on 
the distribution of the Yang-Lee zeros and compare it 
with the spurious zeros of the Skellam partition function, 
originated from the truncation. 

In the next section, we briefly summarize the gen¬ 
eral relation among the probability distribution, partition 
functions, and Yang-Lee zeros. A chiral random matrix 
model and its Yang-Lee zeros are presented in Sec. IIIII 
We demonstrate differences of the truncation effects on 
the Yang-Lee zeros between the models with and without 
phase transition in Sec. El Implications for heavy ion ex¬ 
periments are discussed in Sec. |Vl Section |Vl] is devoted 
to concluding remarks. Detailed expressions for partition 
functions in the chiral random matrix model are given in 
the Appendix. 
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II. GENERAL FRAMEWORK 


We start from experimentally measured data of net- 
baryon number multiplicity distribution T’{N), where N 
is the net-baryon number. In real experiments one mea¬ 
sures the net-proton number AiVp = Np — Np. In princi¬ 
ple one can reconstruct V^N) from 'P{ANp), V{Np) and 
V{Np) [s^. In this study we entirely assume the isospin 
invariance and regard V{ANp) as a proxy of V{N). The 
shape of the distribution depends on the colliding en¬ 
ergies, centrality etc. The net-baryon number can take 
any value as long as it can be packed within the system 
volume. Owing to limited statistics, however, we do not 
observe such states that have too large N far from its 
mean value M. Thus, we define the possible minimum 
and maximum of Wnin and Amax as 

ViN < = 0 

V{N > A„,ax) = 0. (1) 


In thermal equilibrium, probability distribution of the 
net-baryon number in the grand canonical ensemble 
reads, for the fugacity factor A = 


P(T,t/,A,/r) 


Z{T,V,N)X^ 


( 2 ) 


where Z (T, V, fx) is the grand partition function 


Z{T, V, n) = (3) 

and Z (T, V, N) is the canonical partition function 


Z(T,I/,A) = Tr[e-*/^<5^ J. (4) 

Assuming the measured multiplicity distribution is the 
equilibrium one, one finds N dependence of V{N) comes 
from Z(T, V, N)X^. Using the charge-conjugate symme¬ 
try Z{T, V, —N) = Z{T, V,N), one can determine A from 
V{N) and obtain the canonical partition function 

Z{T,V,N)=P{N)X-^. (5) 


Because of limited range of N o, the canonical 
partition function o can be obtained for N G 
[max(0, Amin), Amax] ■ For the energy scan range in RHIC 
experiments, Amin < 0, ie., there are a few events in 
which more anti-protons are observed than protons, ex¬ 
cept for y/SNN = 7.7 GeV where Amin = 1 [3- Note 
that we need only A > 0 thanks to the charge conju¬ 
gate symmetry. Thus, in most cases, one can extract the 
canonical partition function for — Amax A A < Amax- 
This limitation in A also applies to theoretical ap¬ 
proaches such as model studies [1^ |4l|, and lattice 
simulations [13, In the former, the canonical par¬ 
tition functions have been calculated using a projection 
formula 


Z{T,V,N) 


1 

27ri 


dX 


Z{T,V,X) 

XN+l 


( 6 ) 


where integration contour C in complex A plane can be 
arbitrary, but it is convenient to take the unit circle A = 
e*®. Then the formula becomes 

Z{T, V,N) = ^ [ de cos(A6»)Z(r, V, 9) (7) 

27r 

where 9 is related to an imaginary chemical potential 
—ifx/T. In Refs. [1^ EH, 113: thermodynamic potential 
= —pV in Landau theory (4ll | and in chiral quark- 
meson model [13 Ell was used through Z (T, V, 9) = 
g-n(r.e)/T_ Owing to the rapid oscillation in large A, 
it turned out that the numerical integration in double 
precision works up to P{N) ~ 10“^^. 

In lattice QCD simulations, two approaches can pro¬ 
vide the canonical partition functions, i.e., (Bthe fugac¬ 
ity expansion of the fermion determinant |57l | and (ii) 
Hasenfratz and Toussaint method [13. In the fugacity 
expansion, we must diagonalize a matrix whose rank is 
proportional to the lattice spacial volume. This requires 
large computational resource and currently one cannot 
go to simulations on large lattices. In the method (ii), as 
A increases, more accuracy is needed, and consequently 
A cannot go to very large. 

Once the canonical partition function Z (T, V, A) is ob¬ 
tained, one can also construct a truncated grand canoni¬ 
cal partition function as a series in A 


Z^\T, U, A; A„,ax) 


N max 

^ ZipV,N)X^. 


( 8 ) 


Owing to the truncation of the series at A = —Amax 
and Amax; this partition function is an approximation 
of the exact partition function which could be obtained 
if one can take Amax = A* with A* being the num¬ 
ber of net-baryons fulfilling the system volume [ 13 ] • For 
lattice QCD at finite temperature with A^ x At lattice, 
N* = 2A3 [IHQ. Thus, one needs to establish relations of 
physical quantities obtained from the truncated partition 
function ([5]) with those from the exact partition function. 
As seen in the summation running from — Amax; relativis¬ 
tic partition functions contain negative powers of A. The 
suppression of high A contribution to Z cannot be real¬ 
ized by small A. One needs to know large A behavior in 
ZiT,V,N). 

Similar studies on higher order cumulants of the net 
baryon number have been carried out in Ref. (Hj j. in 
which sufficiently large Amax depending on the order of 
the cumulants is shown to be necessary to obtain a cor¬ 
rect value of the cumulants. In this paper, we focus on 
Yang-Lee zeros for the baryon chemical potential. 

The zeros of the partition function in complex chemical 
potential plane can be obtained by solving an equation 


Z{T,V,p) = 0, (9) 


^ In Ref. [ 5 ^ . N* was derived for the quark fugacity series as 
N* = 2NcN^ 
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for complex /r. Owing to the negative powers of the fu- 
gacity, the equation is a polynomial one in A with order 
2N*. For the truncated partition function one needs 
to solve 


= 0 ( 10 ) 


of which the order of the polynomial is 2A^max- 

In the exact case, the roots have both real and imag¬ 
inary part and its distribution in the complex chemical 
potential plane is expected to form a line, which crosses 
the real axis at the transition point in the thermodynamic 
limit. The behavior of the distribution depends on the 
nature of the phase transition. In Ref. , the behavior 
of the Yang-Lee zeros around the CP was studied by us¬ 
ing a chiral random matrix model. The singularity asso¬ 
ciated with the CP appears as a branch point in complex 
fjL plane and its property is shown to be connected with 
the universality. In lattice QCD, the phase transition 
between different Z{Nc) sector in the deconfined phase, 
Roberge-Weiss phase transition, has been recently ana¬ 
lyzed from a view point of Yang-Lee zeros [55l| . In this 
work, we use a chiral random matrix model similar to 
used in Ref. but with an extension to periodic prop¬ 
erty in imaginary chemical potential as it is necessary to 
have integer net baryon number. 

The partition function is written as a polynomial in 
A = Since a complex root Ai is accompanied with 

its conjugate Ai and the charge conjugate symmetry im¬ 
plies 1/Ai and 1/Ai are also roots, only the roots located 
in the first quadrant of the complex /r plane are inde¬ 
pendent. In practice, it is convenient to use Joukowski 
transformation w = A -I- I/A(= 2 cosh Pfi) and reorganize 
the series in terms of w to reduce the number of roots to 
search for. Using a property of Chebychev polynomial 
r„(cosha::) = cosh(na:), one finds 


A^ -I- ^ = 2 cosh{l3fiN) 

= 2Tjv(w/2). 


2rAr(cosh Pfi) 


( 11 ) 


Then Eq. (I8|) reduces to a series expression containing 
only positive powers. After expanding the Chebychev 
polynomial by Eq. (TATI) . one finds. 


^max 

= Z(r, V,N = 0)+ 

n—1 


X 


[n/2] 

E(-F 




(^ - fc- 1)! n-2fc 

k\{n-2k)\ 


( 12 ) 


This formula could be also useful to compare a relativistic 
system with nonrelativistic ones. The roots of uj space is 
easily converted into those in A and /r plane as 

^ = ±cosh"^| (13) 

Taking both signs, one can finds all the roots in the com¬ 
plex /i and A plane. 


III. CHIRAL RANDOM MATRIX MODEL 

In this section, we introduce a chiral random matrix 
model which is an effective model for the spontaneous 
chiral symmetry breaking in QCD. Since this model is an¬ 
alytically solvable in the chiral and thermodynamic limit 
|59l| and analytic expression for the partition function in 
finite volume is known we find this model as the 
most suitable one for the present purpose. An apparent 
shortcoming of the model for applying to the net baryon 
number probability distribution is lack of periodicity in 
imaginary chemical potential, which is a consequence of 
U{1)b symmetry. Thus, we first extend the model to ex¬ 
hibit the appropriate periodicity and the phase structure 
in the imaginary baryon chemical potential. 

In QCD, the partition function has a periodicity 27r/iVc 
in the imaginary quark chemical potential 0q = 0/3, thus 
27r in the baryon number. LQCD simulations have shown 
that there is no phase transition in imaginary baryon 
chemical potential at temperatures below chiral crossover 
temperature and thermodynamic quantities smoothly be¬ 
have as ^ cos 9 [60l - [^ . This fact combined with Eq. © 
implies that the phase transition at large baryon number 
density is encoded in higher Fourier coefficients of the 
smoothly oscillating function. 


A. Partition function and thermodynamics 


We start with a partition function of the chiral random 
matrix model with Ns sites given in 

Zrm = J VX exp |^-^TrAAl'^ det^^ {D + m) (14) 


where a denotes the variance of the random matrix X 
which has Ng x Ng dimension and D is the 2Ns x 2Ng 
matrix approximating the Dirac operator. At T = /r = 
TO = 0, tr is the only dimensionful parameter. We use it as 
a unit of mass in the model and put cr = 1 in expressions 
below. 

The Dirac operator takes the form 


/ 0 iX + iC\ 

At + iC 0 ) ' 


(15) 


The matrix C describes the effect of temperature and 
chemical potential. In Ref. [5^, it was chosen as 


Ck = a^T + ^ 

iN c 


(16) 


for one half of eigenvalues and 

Ck = -anT + ^ (17) 

tiy Q 

for the other halfd, with a and b being the dimensionless 


^ Note that is the chemical potential of the baryon number, thus 
fi/Nc stands for that of the quark number. 
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parameters. 

The linear ansatz for the matrix C (Hsi-dni) accounts 
for the fact that there are the two smallest Matsubara 
frequencies zLttT. This model does not have any thermal 
distribution which gives the fugacity factor thus nor 
the periodicity in imaginary chemical potential, since it 
appears as a result of summation over the Matsubara 
frequencies. In order to make the partition function pe¬ 
riodic, we perform a following replacement 

+ = + (18) 

which gives a periodicity 2ttT in /r/ to the partition func¬ 
tion. Compared to the original linear ansatz, this replace 
does not change anything at /r = 0 but alters the phase 
structure at /r > 2T. 

The phase structure of the model is easily evaluated 
by taking Ng —>■ oo limit. Introducing an auxiliary Nf x 
Nf complex matrix field </> and performing the Gaussian 
integration with respect to X, one obtains the partition 
function 

-Zrm = J V(j3eyii>[-Ngn{(l))] (20) 


1.2 
1.1 

1 

o 

0.9 

H 

0.8 
0.7 
0.6 

0 0.5 1 1.5 2 2.5 3 3.5 

H/Tc 

FIG. 1. Phase diagram of the chiral random matrix model 
with periodicity in the imaginary baryonic chemical potential. 

where 

f = naT (26) 


m=0, 2nd order 
m=0,1st order — 
m=0.05,1st order — 
TCP • 
CP • 



where n((/>) stands for the effective potential. Then the 
partition function can be determined by the minimum of 
the potential, which is evaluated at the saddle point (j)Q 
of the integrand: 


d(j) 


4>=4>o 


= 0 , 


( 21 ) 


and 


lim —ln2^RM = — minn(()f)). 

Vs— »-oo Ns 4> 


( 22 ) 


The saddle point is related to the chiral condensate 
through 




1 9 In Zbm 

NfVi dm 


1 Ns 

NfVi a 


2ReTr((>o 


(23) 

(24) 


where the four dimensional volume V 4 corresponds to Ng 
such that Ng represents the typical number of the instan- 
ton (or anti-instanton) in V 4 . For real m, one expects 
is a real matrix proportional to the unit matrix. There¬ 
fore, the saddle point can be obtained by solving (|^ for 
the potential 




-In 


{ 4 > + m f - sinh ^ 

+ -f^ (a sinh^-i^ | 


(25) 




2 b 

anNc 


(27) 


In the chiral limit m = 0, one finds that t/fg = 1 at 
T = ^ = 0 and a second order phase transition occurs at 
T = l/ina) and /r = 0, where (po continuously approaches 
to zero. Thus, pQ can be regarded as an order parameter 
of the chiral phase transition. 

The parameters in the model, cr, a and 6 , are deter¬ 
mined as follows. The only dimensionful parameter cr 
is estimated to be cr ^ lOOMeV through Eq. (l24)l by 
putting (tjj'fp) ^ 2 fm“^ at T = ^ = 0. Since Tc = l/( 7 ra) 
at ^ = 0, putting Tc = 160 MeV yields a = 0.2. The 
remaining parameter h connects the model to the den¬ 
sity scale. With the linear ansatz for C (HlD-dni) , one 
finds the first order phase transition at T = 0 and 
b^./Nc = 0.528. We follow the choice of Ref. [5^ and 
put b = 0.13, corresponding to the first order transition 
point at fic — 1200MeV, though we do not have the same 
phase diagram as Ref. [59l| owing to the implementation 
of the periodicity (Hiii-dni). 

Figure [ 1 ] shows the phase diagram of the modified ran¬ 
dom matrix model (1251) in the chiral limit and in the pres¬ 
ence of a small explicit symmetry breaking, m = 0.05, 
respectively. In the chiral limit, second order line con¬ 
tinues with decreasing temperature down to T > T 3 and 
^ < /C 3 , where T 3 ~ 0.731Tc and ^3 ~ 4.504 is the loca¬ 
tion of the tricritical point (TCP). Below T 3 , there is the 
first order phase transition line. At finite quark mass, 
the second order line is replaced by smooth crossover 
and TCP becomes CP with slightly decreased temper¬ 
ature and increased chemical potential, Tcp = 0.675rc 














6 



0 

FIG. 2. Behavior of order parameter i^o of a periodic chiral 
random matrix model in imaginary chemical potential 6 = 
Hi/T at nn = 0. 


this form also gives a reasonable thermodynamic quanti¬ 
ties at imaginary chemical potential]! Figure [2] displays 
the behavior of the order parameter (j)Q in the imaginary 
baryonic chemical potential 9 = Hi/T. One sees that 
our parameterization 113 gives the correct periodicity 
27r and expected temperature dependence such as larger 
amplitude at higher temperature below Tc [s!, . Ow¬ 

ing to lack of a ^(3) sector such as the Polyakov loop 
background, this model does not exhibit the Roberge- 
Weiss phase transition at high T. 

At finite the partition function can be expressed 
as [ 5 ^ 


Ns /2 

-Z^RM = ^ 

ki ,^2=0 


Ns/ 2 \ (Ns /2 
ki 


{N-ki-k 2 )l 


XiFi(fci + fc2 - Ns-, 1; -m^Ns){-Nsf^)>^^+^^ 

2 ki / 11 \ 2 k 2 


X (i + Asinh^'^ 


^ (■ 
V 


_ Asinh-^y 

2T) 


(28) 


and ^cp = 4.72, respectively. While these structures are 
the same as those in Refs. [5a,l5^, the apparent singular¬ 
ity at T = 0 in the periodic parametrization significantly 
modifies the phase boundary at low temperature. We 
stress that our purpose in this paper is to explore the 
property of partition function zeros rather than determi¬ 
nation of the phase structure. 

With the parameter set for a, b and cr, we hnd that 


where an irrelevant constant factor is ignored and 
iFi{a,b-,x) denotes the confluent hypergeometric func¬ 
tion. One may directly obtain zeros of this partition 
function, but one needs to expand Z in a series of the 
fugacity A to examine effects of tails in the probability 
distribution function. We put the details in the Appendix 
s and write down only the result for the canonical par¬ 
tition function, for d = |fci — /C 2 I, 


Z{T,Ns,N) = ^ - fc2)!iFi(fci + fc2 - Ns-, 1; -m^Ns){-Nsf^A^ 

ki,k2=0 ^ 1 / V 2 


f ki-\-k2 

E 


fci -I- k2 
k?. 


2(2- A2) 


A2 


fcl+fe2—fes 


X 


k 3=0 
S 

4lV [6 - k 3 )\i 2 k 3 )l 


h 

ks-N 

2 


16 ^ ((5 -I- fca - 1)! _|_ ^2 - fca 


2(2-A^) 


1 k\-\-k 2 — A13 — ^4 


A2 


ki 

ki-N 

2 


ki = k2 

ki yf k2 
(29) 


B. Phase boundary and Yang-Lee zeros 

We compute the Yang-Lee zeros for the truncated par¬ 
tition function m with the canonical partition func¬ 
tion of the chiral random matrix model ([2^ . Taking 


^ Note that in Refs. [SS [59| the coefficients in the temperature 
and chemical potential are absorbed into T and fi. While the 
qualitative phase structure does not depend on the parameters 
in the linear ansatz, it does so when one employs the periodic 
parametrization Gil. 


A^max = Ns in Eq. dni), one recovers the exact grand 
partition function (1281) . The computation of the zeros 
requires a special care in numerical digits as cautioned in 
literature [^, [ 63 . We perform the calculations in 50-300 
digits utilizing FMLIB package in fortran 90. 

Figure [3] shows the distribution of the Yang-Lee zero of 
the periodic chiral random matrix model in the complex 
w plane, for m = 0 and at T/Tc = 0.99. The distribution 
of the zeros is symmetric with respect to the horizontal 
axis because the partition function is an even order poly¬ 
nomial of (jj and a root has its complex conjugate. The 
solid line in Fig.[3]stands for Stokes boundary, which can 
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Reco 


FIG. 3. Yang-Lee zeros of the periodic random matrix model 
in complex lo plane for m = 0 at T/Tc = 0.99. Open symbols 
stand for the zeros in different Ns and solid line indicates the 
Stokes boundary. Branch point is denoted by closed circles. 


be regarded as an extension of the phase boundary to a 
complex chemical potential plane. In the thermodynamic 
limit Ng —>■ oo, it satisfies 


Re 


d(jP 



> 0 , 


ReQ{(j) = (j)o,i) = Ren{(j) = ^ 0 . 2 ), 


(30) 

(31) 


where the first condition ensures the well-defined parti¬ 
tion function at the saddle point of integrand in Eq. (l20l) 
and the second condition denotes the continuity of the 
real part of pressure at the boundary [^. (po^i and ^ 0,2 
stand for the two out of five solutions of the gap equation 
dfl / d(j) = 0 and give the minimum of Refl in both sides 
of the boundary, respectively. The density of the zeros 
increases with Ng and turns into the cut which consti¬ 
tutes the Stokes boundary in the thermodynamic limit. 
This is clearly seen in Fig. [31 There are two branch points 
located on the real axis. Since w = A-|-l/A>0for real /r, 
the one at Rew = Wc = 3.4 > 0 corresponds to the second 
order phase transition point in real while the other 
one, Rew = —634.4 is located on the line Im^/T = tt. 
The Stokes boundary exhibits a closed curve, reflecting 
the periodicity in imaginary ^ and existence of the phase 
boundary at real ^ axis and Im/i/T = tt. 

The phase structure can be more intuitively under¬ 
stood by going to complex /r plane. Figure [Hleft dis¬ 
plays the distribution of the same Yang-Lee zeros as in 
Fig.ia but the zeros in Re/r < 0 are omitted since their 
locations are trivial according to the charge conjugate 
symmetry /r —> The branching point on the hor¬ 

izontal axis indicates the second order phase transition 


point. The Stokes boundary extends to both direction in 
imaginary and ends up at the other branch point. Note 
that the branch points at Im^/T = —tt and tt are essen- 
tially the same because of the periodicity. We refer to 
[691170|| for behavior of the order parameter in complex 
/i plane and related topics. The zeros distribute along 
the boundary and becomes more dense for large Ng, but 
distance to the real axis is not so close for these values 
of Ng. The behavior of the density of the zeros is related 
to a property of the thermodynamic potential which can 
be described by an analogy to electrostatics [s^. In this 
case, ReO can be regarded as the electrostatic potential 
on the (Re/r/T, Im^/T) plane and the normal component 
of the electric field E = — V(Ren) to the Stokes bound¬ 
ary has a discontinuity of which amount is proportional 
to the density of the zeros. We confirmed that in this 
model these discontinuities at large Re/i, where the zeros 
are dense, are much larger than those at small /r, follow¬ 
ing the expectation. Although the density of the zeros 
far from the branching point is a model-dependent fea¬ 
ture dependent on the shape of the Stokes boundary, it 
is governed by the universality near the branch point on 
the real axis as pointed out in Ref. [^ . 

Effects of the finite but small quark mass can be seen 
in the right panel of Fig. 0] where the distribution of 
Yang-Lee zeros for m = 0.05 at the same temperature 
is displayed. Owing the explicit chiral symmetry break¬ 
ing, the phase transition becomes a crossover such that 
the branch point on the real axis moves to above. As 
a result, there are two branch points of which are com¬ 
plex conjugate each other. The same thing occurs also 
to the branch point at Im^/T = ±7r. Here we emphasize 
that these complex singularities are, albeit unphysical, 
indicating existence of a chiral phase transition in the 
chiral limit. These are also known as Yang-Lee edge sin¬ 
gularities [13 . The critical point at finite density (See 
Fig.® is realized by coalescence of the branch points 
close to real ^ axis when temperature is decreased [s^ . 
As seen in Fig. |3]-right, the Yang-Lee zeros are fairly on 
the boundary line and exhibit expected behaviors. 


IV. YANG-LEE ZEROS FROM TRUNCATED 
PARTITION FUNCTIONS 

As described in Sec. im the connection of net baryon 
number multiplicity distribution m with the recon¬ 
structed grand partition function ([5]) could potentially 
enables us to extract the Yang-Lee zeros from experi¬ 
mental data. Since the results presented in the previous 
section correspond to Vmax = N*, i.e., no information 
on the exact partition function is lost, we need to evalu¬ 
ate whether one can obtain the correct distribution of the 
Yang-Lee zeros when the fugacity expansion is truncated. 
Furthermore, even if one starts from a partition function 
which does not exhibit any phase transition such as an 
ideal Boltzmann gas, the truncation produces the zeros 
of partition function because it is a polynomial of order 







Re|X/T 


Re|X/T 


FIG. 4. Yang-Lee zeros of the periodic random matrix model in complex /r plane. The left panel corresponds to the case of 
Fig. [3] Right panel shows the case with a finite but small quark mass, m = 0.05 at the same temperature. 


Nmax- In this section we investigate in detail the effects of 
the truncation on the distribution of the Yang-Lee zeros. 

A. Random matrix model 

Figures [5] and [6] display the distribution of the Yang- 
Lee zeros from the truncated partition function of the pe¬ 
riodic chiral random matrix model for various fVmax and 
m = 0. 050 Hereafter we set Ng = 60. We confirmed the 
following results does not depend on the specific choice of 
Ng. We plot only the first quadrant in complex ^ plane 
according to the symmetry structure of the distribution. 

The left panel in Fig. [S] shows the case of T = Tc, at 
which transition is of crossover type as seen in the branch 
point at (Re/i/T, Im^/r) ~ (3,7r/4). For N^ax = 60 = 
Ng, the zeros are located on the Stokes boundary (dashed 
line). Reducing A^max by one, i.e., removing Z{N = 60) 
from the series, one sees a drastic change in the distri¬ 
bution. The distribution of the zeros at large Re^ and 
Im^ splits into the two lines, but the rest of the zeros 
remains unchanged. Further reduction of Wnax substan¬ 
tially modifies the distribution such that the splitting 
occurs closer to the edge closer to the real /r axis. Nev¬ 
ertheless, up to Nmax = 21, the edge of the distribution 
which is the closest Yang-Lee zero to the real p, axis re¬ 
mains the same. Beyond iVmax = 20, the distribution no 
longer holds the information on the exact Yang-Lee zeros 
thus the apparent relation to the phase boundary is lost. 


^ Note that Tc is defined for m = 0. Thus it is slightly lower than 
the chiral crossover temperature for m = 0.05. 


The behavior with respect to changing N^ax does not 
depend on temperature or corresponding phase transi¬ 
tion. In the right panel of Fig. [SJ we plot the result 
of the same analysis for T = Tcp = Q.&TbTc where the 
branch point appears on the real axis, indicating the crit¬ 
ical point. Reflecting the location of the branch point, 
the edge of the distribution also become closer to the 
real axis compared to the crossover case. The edge is 
stable against decreasing Wiax down to N^ax = 19, 
then it starts to deviate slowly when decreased further. 
This is so also in the case of first order phase transition 
(T = 0.6Tc) depicted in Fig. [6l The branch point is hid¬ 
den in unphysical Riemann sheets [t^ and the edge is 
very close to the real axis. 

We also note that there is always a zero at Im^/T = 
TT when A^max is odd. These zeros look special since it 
corresponds to negative real axis in both complex A and 
uj plane. However, this is a mathematical consequence 
because in this case the truncated partition function (11211 
is an odd order polynomial, thus it has at least one real 
root. As seen in Figs. [5] and [SI it becomes the edge of one 
of the lines bifurcating from the exact Yang-Lee zeros. 

These results indicate the stability of the edge does 
not depend on the detail of the phase structure, although 
the location of the edge seems to be connected with the 
shape of the Stokes boundary which is model dependent 
through the /i dependence of the partition function. In 
particular, the present results are obtained by employing 
the periodicity (1191) in the random matrix model which 
does not correctly take into account degrees of freedom 
with baryon charges . We note that this modification 
causes unphysical behavior in thermodynamics, such as 
negative Z{T,V,N) at some small N at low T, which 
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FIG. 5. Distribution of the Yang-Lee zeros from the truncated partition function of the periodic chiral random matrix model 
for m = 0.05. Left and right panels stand for T = Tc and T — Tcp, respectively. 
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FIG. 6. Same as Fig. [5] but for T = 0.6Tc 


FIG. 7. Distribution of Yang-Lee zeros for a lattice QCD data 


presumably reflect the unusual curvature of the phase 
boundary in Fig. [T] Therefore, we note that the shape 
of the distribution itself might not be relevant for realis¬ 
tic situations. Nevertheless, below we shall see that the 
stability of the edge is specific to the case with a phase 
transition. 


B. Lattice QCD 

Figure [7] displays distribution of the Yang-Lee zeros 
above Tc calculated in lattice QCD simulation via the 
canonical method . While calculations in the confine¬ 
ment phase is still numerically difficult thus we do not 
see clear indications of a phase transition at low T, the 
Roberge-Weiss (RW) transition provides us a well- 
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defined phase transition in high temperature quark-gluon 
plasma phase, though at imaginary chemical potential. 
In this figure, the data are calculated on 8^ x 4 lattices 
and /? = 1.89 which corresponds to T/Tc ~ 1.94. A more 
detailed analysis in lattice QCD with different lattice se¬ 
tups can be found in Ref. [55|. Since quark mass is heavy, 
the calculation is not relevant for chiral phase transition. 
The RW transition is regarded as a transition from one 
Z{3) sector to another one when single quarks can be 
excited owing to deconfinement and is known to exhibit 
a first order phase transition at Im^Xg/T = ±7r/3. In 
terms of baryon chemical potential, the transition lines 
reduce to Im/i/T = ±7r which is shown as a dotted line in 
Fig. [71 A brief explanation of the Roberge-Weiss phase 
boundary can be found in Appendix |B1 Since it is hard 
to compute Z{N) near N = N*, the canonical approach 
in lattice QCD lacks large N contribution when one con¬ 
structs the truncated partition function (ITTl) . One sees 
that in Fig.jTjthe behavior of the distribution of the Yang- 
Lee zeros against changing A^max is similar to that of the 
random matrix model, despite the completely different 
origin of the phase transition. Therefore, we expect the 
similar splitti^ behavior of the distribution also appear¬ 
ing in Refs. is also due to the truncation effect. 

Indeed, Ng and fVmax dependence of the Yang-Lee zero 
shown in Ref. agrees with the truncation effects dis¬ 
cussed here. We expect that the bifurcation of the zero 
starts at large Re/i by improving the fugacity expansion, 
but one needs to take iVmax = N* to completely produce 
the Yang-Lee zero along the transition line. In the RW 
transition where the transition point at luifi/T = ±7r, the 
edge of the distribution is the closest zero to the imagi¬ 
nary axis. One sees that this point is also stable against 
changing A^max- This fact suggests that the stability of 
the edge is not specific to the random matrix model or 
chiral phase transition but might be a general property 
of the distribution when partition function is truncated. 


C. Skellam distribution 

Finally we examine a model without phase transition 
in order to check whether the stability of the edge is spe¬ 
cific to phase transition or not. We employ the Skellam 
distribution 0 of which probability distribution of the 
net baryon number N is given by 

/AT \ 

P^N) = lN{2A/N^)e-^^‘^+^s) ( 32 ) 

where Nb and Nb denote the thermal averages of the 
numbers of baryons and anti-baryons, respectively. The 
mean M and variance of the distribution are given 
by M = Nb — Ng and = Ng + Ng, respectively. 
For Ng = Ng, The distribution becomes symmetric and 
the argument of the modified Bessel function In{x) is 
reduced to 2Nb = This distribution can be derived 
from non-interacting Boltzmann gas [l^ , thus the canon¬ 


ical and grand canonical partition functions read 


ZiN) = 

Z(A) = exp ^ 



(33) 

(34) 


where the temperature and volume dependence is en¬ 
coded in (T^. Obviously the grand partition function (l34l) 
does not have any roots thus no phase transition exists. 
When one constructs the truncated grand partition func¬ 
tion ([5]) from the canonical partition function (1331) , how¬ 
ever, there exist complex roots. Consequently, one might 
see these spurious zeros even if the system does not have 
any phase transition, when one constructs the partition 
function through the fugacity expansion. 

Here we investigate such spurious zeros from the Skel¬ 
lam partition function (I33|) such that it has the same 
variance with the random matrix model at Ng = 60, 
T = Tc and m = 0.05 of which distribution of the Yang- 
Lee zeros is displayed in Fig. |SJ Since the information on 
the phase transition is encoded in the tail of the prob¬ 
ability distribution P{N), the Skellam distribution with 
the same variance serves a useful reference distribution 
[^ . The probability distribution of the random matrix 
model and corresponding Skellam distribution are shown 
in Fig. |5]-left. Both distributions almost agree for small 
N, according to the same tr^, but the deviation appears 
in the tail of the distribution with tiny probability. 

Figure [Sl-right shows the distribution of zeros of the 
truncated partition function for the Skellam distribution 
with cr^ = 0.365. Except for a splitting of the distribution 
for A^max >70 which is similar to those in the random 
matrix model, the distributions consist almost parallel 
lines moving to large real /r direction as A^max increases. 
This behavior reflects the fact that all the zeros go away 
to infinity as A^max oo since the exact grand partition 
function does not have roots. Remarkably, the edges of 
the distributions also move together with the rest of ze¬ 
ros, in contrast to the random matrix model and lattice 
QCD. Furthermore, one notes that the distributions for 
A^max < 20 in the random matrix model, shown in Fig.jSJ 
resemble those from the Skellam distribution. This ob¬ 
servation indicates that the stability of the edge against 
Nmax is a consequence of the existence of phase transi¬ 
tion and information on the phase transition is lost for 
too small A^max- 


V. DISCUSSION 

A. Comparison with Amax for cumulants 

In the previous section, we have shown that the edge 
of the distribution of the Yang-Lee zeros remains un¬ 
changed when the tail part of the canonical partition 
function is missing. In practice, this property gives im¬ 
plications for necessary statistics in heavy ion studies of 
the net-baryon number fluctuations and in lattice QCD 
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FIG. 8. Left: P{N) for the random matrix model at T — Tc and m = 0.05, and the corresponding Skellam distribution. Right: 
Distribution of Lee-Yang zeros for a truncated Skellam partition function. 


TABLE 1. Ymax necessary for reconstructing i-th order cu- 
mulants and edge of the Yang-Lee zeros from Z{N) in the 
random matrix model at T = Tc and m = 0.05. 


Ns 

■1 ’max 

’ max 

’ max 

^max 

60 

3 

4 

6 

21 

80 

3 

5 

6 

26 

100 

4 

5 

7 

30 


calculations. Since the sufficient /Vmax to see the stable 
edge depends on the system volume, here we compare 
it with corresponding for f-th order cumulants. 

Here we consider only even order ones for net-baryon 
number at /r = 0, since we are looking at Z{T,V,N) 
rather than P{N) which becomes asymmetric with re¬ 
spect to N at fi > 0. Thus the first central moment 
SN = N — (N) = N. The second, fourth and sixth order 
cumulants c„ (n = 2,4, 6) read 

C2 = mf) (35) 

C4 = {{6N)^) - 3{{SNff (36) 

eg = {{6Nf) - 15{{dN)^){{dNf) + 30{{6N)Y (37) 


The property of the higher order cumulants of net- 
baryon number probability distribution for changing 
Nmax was studied in Ref. by using a chiral quark- 
meson model. For sufficiently large volume, it was shown 
that Nmax for the cumulants approximately scale with 

w. 


We summarize the values of each A^max in Table ID 
The calculations are done for T = and m = 0.05 in 
the random matrix model. Owing to the narrow Z{N), 


even the sixth order cumulant for Ng = 100 only re¬ 
quires A^mix = 7,i.e, —7<A^<7to reconstruct it from 
Z{T,V, N), while the edge of Yang-Lee zeros demands 
Nmax = 21. The small A^max implies the system vol¬ 
ume is not large enough to exhibit Vv scaling regime of 
the cumulants. This can be understood from the small 
value of in the random matrix model calculations. The 
rapid decay of P{N) give a rather weak dependence of 
Nmax for the higher order cumulants. For a sufficiently 
large volume, one expects that P{N) resembles Gaussian 
near the peak, while the probability distribution (Fig. [8]) 
has a sharp peak. Thus we cannot assess the value of 
Nmax needed in a realistic situation relevant for heavy 
ion collisions. Moreover, the baryon number carried in 
this model is not a physical one, as mentioned above. All 
we can say is that one may need much more statistic than 
higher order cumulants. 


B. Skellam distribution for large volume 


In the fluctuation measurements at RHIC, observed 
P{N) can be well described by the Skellam distribution 
and deviation from the Skellam distribution exists in the 
tail, resulting in higher order cumulants different from 
the Skellam case. The obtained variance reaches tr^ ^ 10 
[H, [zlj at the most central bin. Thus it is instructive to 
give a reference for the distribution of the spurious Yang- 
Lee zeros based on the Skellam distribution. Here we pick 
up the data for y/s^N = 7.7 GeV at the most central bin, 
which gives cr^ = 14.89 and M = 14.41 with available bin 
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from N = 0 to N = 340As mentioned in Sec. im one can 
construct Z{N) from — 34 < Af < 34 according to charge 
conjugation symmetry (4^ . In the Skellam distribution 
for (T^ = 14.89, we find that iVmlx = 13,20 and 26 for 
second, fourth and sixth order cumulant, respectively. 
Note that these fVmax apply to the cumulants at ^ = 0. 
The data does not have enough statistic for the sixth 
order cumulant at freeze-out /i. 

We plot the distribution of the Yang-Lee zeros for the 
constructed Skellam distribution with = 14.89. The 
basic feature is the same as the small case (Fig. [8|). 
The line of zeros moves toward infinity as iVmax increases. 
One sees that some zeros below iV = 30 appear on the 
imaginary ^ axis, which corresponding to the unit circle 
in complex fugacity plane. This means that, for a large 
volume case, the zeros can appear on the imaginary axis 
when the tail of the Z{N) is not provided. In the Skellam 
distribution, these zeros can be obtained directly by look¬ 
ing at the truncated partition function on the imaginary 
^ axis, for A = e*®. 


VI. CONCLUDING REMARKS 

In this paper, we present analyses on partition func¬ 
tion zeros which can be obtained from a truncated series 
of the fugacity expansion. By solving an extended chi¬ 
ral random matrix model which has a periodicity in the 
imaginary chemical potential, we compare the exact lo¬ 
cation of the Yang-Lee zeros and those obtained from 
the truncated series. We found that the edge of the dis¬ 
tribution of the zeros is insensitive to the truncation of 
higher order terms in the fugacity expansion to some de¬ 
gree. We found the similar behavior in lattice QCD at 
high temperature in the context of the Roberge-Weiss 
phase transition. This observation indicates that those 
higher order terms may have limited influences in search 
for the location of the phase boundary in lattice QCD 
calculations and heavy ion experiments. Although the 
distribution of zeros exist in systems without phase tran¬ 
sition, due to the truncation, the zeros closest to the real 

axis are stable against truncation if the system has a 
phase transition or crossover. The spurious zeros in the 
Skellam distribution moves toward infinity against the 
truncation. Therefore, one can distinguish whether the 
distribution is related to the phase transition or not by 
looking at the stability of the edge of the distribution 
against the truncation. 

Although the information on the Stokes boundary is 
lost in the case of too small A^max, we expect that it 
does not mean that all the relevant information on the 
phase transition gets lost in the truncated partition func¬ 
tion. This expectation follows from the fact the sixth and 
higher order cumulants at ^ = 0 should be influenced by 
the phase transition and the truncated series is still able 
to reproduce them. It would be interesting to see how 
the distributions of the zeros in the small Vmax cases in 
Figs. [SHni are related to the remnant of the phase transi¬ 
tion. 

The order of the truncation in the fugacity series to ob¬ 
tain the stable edge of the Yang-Lee zeros, Vmax, is nev¬ 
ertheless found to be much larger than those for higher 
order cumulants. We cannot make a quantitative assess¬ 
ment on realistic values for heavy ion experiments due to 
the lack of connection in the model to the real world. We 
hope that such an estimate becomes feasible in the near 
future. 


^max 

^rkenam(0)=2 ^ /at ) COS iV0, (38) 

.^max 

which converges into Eq. (1341) with oscillations giving ze¬ 
ros on imaginary /i. 


® The data for A = 0 and A = 34 have only 1 event. 
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Appendix A: Derivation of the canonical partition 
function in a chiral random matrix model 


Appendix B: Roberge-Weiss transition as a thermal 

cut 

In this appendix, we give a brief explanation of the cut 
arising from the Fermi distribution function and apply it 
to the Roberge-Weiss transition in QCD. 


1. Thermal cut in free Fermi gas 


In the following, we derive an analytic expression for 
Z[T, Ns, N) from Eq. (05)) . First we rewrite the /j, depen¬ 
dent part in terms of the fugacity A = . 

Since 


Asinh^+x)"''^ 
2T J 

{A sinh — i 

V 2T 

2k2' 

(Al) 

(A^sinh^ATl^ 

ki+k2 



cos[2(fci - 

- k2)4‘] 

(A2) 


The thermodynamic potential of the free Fermi gas is 
given by 

j ln[l + + (m ^ -m) (B1) 

where Ep = When the chemical potential /i 

has an imaginary part, /r/ = 9T, the imagary part gives 
a phase in front of the Boltzmann factor: 


where 


I _|_ g-PiEp-IJ.) = 1 _|_ g*Sg-/3(Bp-/^K) 


(B2) 


tan(/) = (^Asinh^) , (A3) 

and imaginary part vanishes after summation over ki and 
^ 2 , using the Chebychev polynomial 

Tk^-k 2 (cos 2(j)) = cos[2(fci - k 2 )(j)] (A4) 

and 

cos 2 (j) = sinh ^ “ l) / sinh ^ + l) . (A5) 


where /i = -I- i/x/. Therefore, for 9 = ±7r, the phase 

gives —1 and the thermodynamic potential has a loga¬ 
rithmic cut at 9 = ±7r and m < fiR < oo. The anti¬ 
particle term also gives the cut symmetric with respect 
to the imaginary axis. In Ref. [t^ . it is pointed out that 
the branch point singularity limits the convergence ra¬ 
dius when one tries to analytically continue the results 
in the imaginary chemical potential to the real one. Since 
this cut originates from the Fermi distribution, the same 
analytic structure appears in chiral models with fermions 


we have the partition function Z as 


2. Roberge-Weiss transition 


NJ2 

E 

ki ^k 2—0 


Ns/2\ (Ns/2 
fci 


(Ns-k^-k 2 )\ 


XiEi(fci + fc2 - Ns, 1; -Tr?Ns){-Nsf'^A^I^f^^'^^ 


X 




2(2-A2) 


ki+k 2 


/A + A-1-2(A2 + 2)/A2\ 

U + A-i-2(A2-2)/A2j ■ 


(A6) 


Expanding the Chebychev polynomial by the following 
expression 


In QCD at high temperature, quarks are deconfined 
and have a light mass owing to chiral restoration. Since 
the deconfinement can be expressed as a breaking of 
Z{Nc) symmetry, it is useful to resort to chiral effec¬ 
tive models with the Polyakov loop background 
which successfully describe the Roberge-Weiss transition 
[S EE ■ Then, the relevant leading single quark con¬ 
tribution to the thermodynamic potential reads 

O _ ^ ^ ^ lull -I- — 

y (2^)3^^^^ + ^^" J (B3) 

+ (/Xg —>■ —Mgi d) —>• $) 


1 


n = 0 


Tn{x) 




fe=0 


{n + k — 1)! 
{n-k)\{2ky. 


{l-xf 


n > 1, 


(A7) 

and using binomial expansion in the third line of (IA6I) , we 
can express Z in terms of A + A“^. We obtain Eq. 
applying the projection (jS]). Note that maximum power 
of A is given by Ns- 


where fig = fi/3 is the quark chemical potential and $ 
is the the expectation value of the Polyakov loop. For 
antiquark contribution, the conjugate <1’ couples to the 
thermal distribution. At the imaginary chemical poten¬ 
tial, Polyakov loop $ acquires a complex phase One 
may express $ = Then the thermodynamic con¬ 

tribution becomes 










14 


The phase of the Polyakov loop (p varies as a func¬ 
tion of the imaginary quark chemical potential 9q. The 
Roberge-Weiss transition at 6 q = tt/S can be understood 
as a transition from one 2^(3) sector with (p = 0 to an¬ 
other one tp = —‘I'KjZ [13. Then, the coupling between 
p and dq gives the prefactor —1 in front of the Boltz¬ 
mann factor. Moreover, |$| ^ 1 in the deconfined phase 


and the prefactor 3 allow this function to have the sin¬ 
gularity at Pq^R = 0. This feature gives the cut drawn as 
RW transition line in Fig. [71 A derivation based on the 
Gaussian P{N) can be found in Ref. [1^. In the confine¬ 
ment phase where |$| ^ 0, this term is suppressed and 
the thermal cut from the quark does not appear. 
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